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Fractional Order Modeling of Atmospheric Turbulence - 
A More Accurate Modeling Methodology for Aero Vehicles 



Abstract: 


This presentation has been prepared for the Aerospaee Controls and Guidanee Systems Committee 
meeting, whieh takes plaee in Cleveland Ohio on Oet. 14-17, 2014. The presentation eovers a 
reeently developed methodology to model atmospherie turbulenee as disturbanees for aero vehiele 
gust loads and for eontrols development like flutter and inlet shoek position. The approaeh models 
atmospherie turbulenee in their natural fraetional order form, whieh provides for more aeeuraey 
eompared to traditional methods like the Dryden model, espeeially for high speed vehiele. The 
presentation provides a historieal baekground on atmospherie turbulenee modeling and the 
approaehes utilized for air vehieles. This is followed by the motivation and the methodology utilized 
to develop the atmospherie turbulenee fraetional order modeling approaeh. Some examples eovering 
the applieation of this method are also provided, followed by eoneluding remarks. 
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Background 

• Traditionally atmospheric disturbance applied to air vehicles structural 
component designs to handle loads induced by severe gusts. 

• Collection of atmospheric turbulence data: 

— Gust accelerometer data collected by NASA & Air Force (1956-75) to 
establish a statistical data base (Ehemberger and Love, 1975) 

— Gust velocity & directionality measurements by NASA in late 80 ’s 

— Gradients of temperature transients over 4 years in early 90 ’s by Concord. 

• For propulsion system the concern is freestream flow velocity, temperature, and 
pressure disturbances entering the inlet and angularity of the flow (coupled w/ 
additional pitch, yaw and roll disturbances), which contribute to inlet unstart. 

• Today for the High Speed Project (HSP) the concern is gust loads and 
AeroPropulsoServoElasticity (APSE) - coupling of propulsion and structural 
dynamics. 
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Background 

Historically conditional probability analysis is applied; that event y will take 
place provided that a certain magnitude atmospheric turbulence E is encountered 
(based on atmospheric data) 

P(_y\E) = P(E)xP(y\E) 


Exceedance rates are typically developed for aircraft designs based on 
probabilistic models (specific to a particular design) 
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Soreide & Tank Validation of a Statistical Atmospheric Model Using 747 
Vertical Gust Exceedance Datas, AIAA, 1997 
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Background 



Atmospheric Turbulence models are based on the Kolmogorov spectrum - 
spectrum originally developed by Tatarski (1961). This spectrum has an energy 
that approaches infinity at low frequencies - Kolmogorov; Russian scientist first 
analyzed atm. turb. in 1940’s . „ 





Kolmogorov 
Spectrum 1961' 


von Karman spectrum, 
1996 by Soreide & Tank 
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Wind and Potential Temperature Spectra as 
reported by Nastrom (1985). For clarity, meridional 
wind and temperature spectra have been shifted one 
and two decades to the right, respectively. 
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Frequency, Hz 

Acoustic wave velocity spectral comparisons 
for the Kolmogorov and von Karman spectral 
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Background 

• The slope of the atmospherie turbulenee speetrum ehanges at about 400 km 
turbulenee length. 

• Speetrum density of turbulenee inereases with deereasing wave number 
(increasing length of turbulence). 

• Spectral also increases with increased turbulence levels - increased eddy 
dissipation rate, 8, (energy/mass/time). 



6 = v^v^ / 


V - kinematic viscosity 
v[ - velocity fluctuation in atmospheric 
region of size / 




"k 


8 specifies the rate at which energy in 
the turbulence breaks down when a 
critical Reynolds number is exceeded 
until it completely dissipates, in the 
process leading to increase wind gusts. 

The assumption is that in worst case 
turbulence a critical Reynolds no. is 

L= 0.07size of turbulence 
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Background 



The Soreide and Tank von Karman spectrum can be used as the amplitude of sinusoids 
at different frequencies to develop a time domain spectrum of certain atmospheric 
turbulence intensity of eddy dissipation rate, s, and integral length scale, L. 
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Alternatively, the Dryden model was developed to 
simulate atmospheric turbulence (Hoblit, 1988) 
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The Dryden model is 2nd order instead of 5/3, 
which underestimates turbulence by 7 dB per 
frequency decade. Exact fractional TF not possible mathematically. 

“ Underestimation acceptable for lower speed vehicles, but can be problematic for 
higher vehicle speeds. 
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Background 

Comparison between von Karman model (exaet) and Dryden model (approximate) 
— Underestimation of speetral eomponents, 7 dB/dee deviation after 3db natural 
frequeney. 
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Fractional Order Modeling Motivation 

• Because the Soreide - Tank model is more difficult to use and because the 
Dryden model significantly underestimates atmospheric turbulence for supersonic 
vehicles the need was identified to develop an accurate fractional order 
turbulence model in frequency domain. 

— This model acts as a low pass filter and only unity amplitude sinusoids are 
needed as inputs to simulate the atmospheric turbulence spectrum associated 
with turbulence. 

• Also, instead of probabilistic models like exceedance rates, the motivation is to 
simulate certain worst case atmospheric turbulence and understand how the 
vehicle model or the controls will perform. 

— The reason is based on the assumption that if a vehicle encounters certain 
atmospheric turbulence conditions, breakdown due to an eddy dissipation rate 
associated with that turbulence (that generates peak gusts) disturbances will be 
a certainty. 

— The new model can still be used in conjunction with probabilistic calculations 
like exceedance rates. 
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Fractional Order Modeling Approach 





Estimate a fractional order TF utilizing order pole - zero pairs (for fractional order 


less than 1) - ex. TF for 1 pole-zero pair, G 


(s/0)Z +1) 
(s/ 0 )p+l) 



• The idea is as the density of pole zero pairs increase the stair case estimate will 
approach a straight line. 


• The challenge is to properly locate the frequencies of the poles and zeroes and the 
intersections of these lines so that the stair case estimate is symmetrically located on 
top of the fractional order slope 

(e.g., dl— d2...dn ) ^(®pl”®pi-l> ®zl”®zi-l’ ®Hpi, ®Hzi) 
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Fractional Order Modeling Approach 

1. Change units of spectrum to more convenient units (i.e., from (m/sec) V(cycle/m) to 
(m/sec)/Hz and relate the wavenumber to the speed of the vehicle. 


Kolmogorov 
Spectrum 
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Spectrum 
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Fractional Order Modeling Approach 

2. Develop a fractional order circuit model for the each type of disturbance, t, based on 
the Von Karman spectral and derive the circuit parameters of the model. 
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Circuit Equivalent of Atmospheric 
Disturbances 


Resistance: 




Capacitance: 


C = 


1 




Natural Frequency: co^ = 
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con 




Gain: 


K, =f{s,L) 


v=5/3, r=l/3 for wind gusts and I /2 for Temp. & Pres., q=xr 

Type of disturbance t - wind gust longitudinal, wind gust transverse, temp., pres. 
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Fractional Order Modeling Approach 



3. Verify that the von Karman type circuit model closely approximates the 
Kolmogorov model, apply scaling as necessary for proportional gain. 



Circuit approximation works sufficiently well. However, its filter dynamics cannot be 
simulated because the circuit is fractional order. 
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Fractional Order Modeling Approach 

4. Instead circuit model and its parameters are used as the basis to derive formulations 
for fractional order TF fit based on first order pole-zero pairs. 
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Derived TF of fractional order approximation: 
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Calculation of other supporting parameters covered in literature: 

1. Atmospheric Turbulence Modeling for Aero Vehicles; Fractional Order Fits, NASATM - 2010-216961 

2. Modeling of Atmospheric Turbulence as Disturbances for Control Design and Evaluation of High Speed 
Propulsion Systems. Journal of Dynamic Systems, Measurement, and Control 134:2, 021009, Jan-2012. 


National Aeronautics and Space Administration 

Fractional Order Modeling Approach 

Example TF fits based on the derived model of first order pole-zero pairs, with 4 poles 
and 3 zeroes (1 pole followed by 3 pole - zero pairs): 




Figure C.9. — Von Karman acoustic wave velocity spectral due 
to temperature gust and its TF fits for different integral scale 
lengths (s = 8.6e-5 m^/sec^). 



Figure C.10. — ^Von Karman acoustic wave velocity spectral due 
to temperature gust and its TF fits for different eddy 
dissipation rates (L=762 m). 


Note that differences in the integral scale lengths only impacts the low frequency 
range, where typical control system designs should have no problem to attenuate 
disturbances in this range. 
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Fractional Order Modeling Approach 

With zoom -in: The staircase like approximation can be seen 
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Fractional Order Modeling Approach 

For control systems design this fractional order filter model is simplified by 
eliminating the integral scale length (using L=762 as the basis - considered standard 
in airplane industry according to Tank) 


Wind gust longitudinal 

C (-W-70-2/9 (s/9.2 + 1)(5/55.0 + 1)(s/335.5 + 1) 

^ (5/1.46 + 1)(5/30.1 + 1)(5/85.7 + 1)(5/1593.1 + 1) 

Temperature 

C -9213^2/ 6 L/ 33.0 + 1X5/45.6 + !)(../ 602.4 + 1) 

(5/1.1 + 1)(5/25.1 + 1)(5/109.8 + 1)(5/816.3 + 1) 

Longitudinal gusts due 
to temperature 

^ ^^0.2/6 (M-1) MyR (^/33.0 + 1)(^/45.6 + 1)(5/602.4 + 1) 

^ j ao (5/1.1 + 1)(5/25.1 + 1)(5/109.8 + 1)(5/816.3 + 1) 


• Similar for transverse wind gusts and pressure TFs 

1 . Modeling of Atmospheric Turbulence as Disturbances for Control Design and Evaluation of High Speed 
Propulsion Systems. Journal of Dynamic Systems, Measurement, and Control 134:2, 021009, Jan-2012. 
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Considerations for Atmospheric Turbulence Specifications 

• Good Source of information in this area for a statistical approach is the publication 

“Atmospheric Disturbance Environmental Definition,” by Tank W G., NASA CR-195315, 
1994. 

• This area also covered (worst case - non statistical) in article “Modeling of 
Atmospheric Turbulence as Disturbance for Control Design and Evaluation of High 
Speed Propulsion Systems” ASME Journal of Dynamic Systems, measurement and control, 
2011. 

— Assumes that if a vehicle encounters an atmospheric turbulence, a critical Reynolds 
number will be exceeded and an eddy dissipation rate associated with the severity 
of the turbulence and corresponding wind gusts will be experienced. 

— This seemingly relaxes the necessity for statistical exceedance rates, which are 
based on pure chance and doesn’t take into account the weather conditions a 
vehicle is flying in at any given moment in time. 
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Considerations for Atmospheric Turbulence Specifications 

Continued from previous pg. 

“ Assumes worst case turbulence is a combination of wind gusts, temperature and 
pressure disturbances. Even tough, this is left to the user to define what is worst 
case. 

“ Spectrum frequency range {/range — depends on the smallest length scales 

5 

of turbulence (25m), before viscous dissipation starts to become significant. 

— Spectrum sinusoid density chosen such maximum wind gust will not exceed a 
maximum (like 100 miles/hour - but user defined). 

— Spectrum sinusoids equally spaced, or purposely coincide with known vehicle 
modes (guidelines provided, but user defined) 
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Applications 

• Applied to the control design of a mixed compression inlet 



Atmospheric dist 



(j 100(5 /27r5 + l)(5/ 27118.8 + 1) 

” 5(5 / 27114.5 + 1)(5 / 271104 + 1) 

(52 /(27t145)2 + 0.985 / 27t145 + 1) 
(52 /(27t5000)2 +1.45/ 2715000 + 1)) 
(52 /(27t250)2 +1.35/ 271250 + 1) 
(52 /(27t 5500)2 +1.45/ 2715500 + 1 
(52 /(27t350)2 + 1 .365 / 271350 + 1) 
(52 /(27 t 6000)2 + 1 .45 / 2716000 + 1) 
(52 /(27t450)2 + 1 .05 / 271450 + 1) 
(52 /(27 t6500)2 + 1 .45 / 2716500 + 1) 
(52 /(27t550)2 + 1.05/27t550 + 1) 
(52 /(27t7000)2 + 1 .45 / 2717000 + 1) 
(52 /(27t650)2 + 0.605 / 271650 + 1) 
(52 /(27t7500)2 + 1 .45 / 2717500 + 1) 


Figure 12. — Feedback control diagram of inlet shock position system via bypass door actuation. 
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Applications 

• Applied to AeroPropulsoServoElasticity (APSE) 
propulsion and structural vehicle dynamics. 



closed loop coupling^^f 



Propulsion System Thrust 





/ 





Closed-loop diagram of APSE simulation 


APSE wing velocity i- 
ASE wing velocity 
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Wing velocity at vicinity of engine inlet for 
APSE system and for ASE alone 
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Applications 





Applied to assessments of vehicle thrust forces spectrum for different eddy 
dissipation rates - engine speed control, but no inlet shock control (external 
compression inlet) 
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• Covered a new modeling approach for atmospheric turbulence for aero vehicles. 

“ Approach models turbulence in their natural fractional order form. 

— Approach is more accurate than traditional models like the Dryden model, which is 
widely used by industry and military - especially for high speed vehicles. 

— Probabilistic exceedance analysis can be conducted with this new approach as in 
the past. Even though, this model was developed for deterministic analysis in mind 
(i.e. for vehicle designed to handle specified worst case turbulence conditions). 

— Considerations for atmospheric turbulence specifications also provided in 
this approach. Even though, this is done independently of the methodology. 

— US patent has been filed. 


Conclusions 
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Thank You! 


I I 



Questions? 
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